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ABSTRACT ol O ^ V ° 

When the forces resulting from viscosity and heat injection are described 
in terms of the thermal 4-potential introduced in the preceding paper, it is pos- 
sible to derive a scalar equation of motion that has the form of a generalized 
Hamilton-Jacobi equation. This equation involves not only the electromagnetic 
and gravitational fields, as in the single-particle case, but also the thermody- 
namic properties of the fluid as characterized by the specific enthalpy, whose 
role is analogous to that of the gravitational potential, and by the thermal 4- 
potential, whose dynamical effects are analogous to those of the electromagnetic 
4-potential. 

The formalism employs a generalization of the canonical particle momentum 
that includes the thermal 4-potential as well as the electromagnetic 4-potential. 
This generalized canonical momentum can be represented in terms of three 
scalar functions by means of a generalized Clebsch Transformation. One of 
these scalar functions is Hamilton's Characteristic Function, which is the un- 
known in the generalized Hamilton-Jacobi equation. The other two functions, 
which are called the vorticity invariants, determine the intrinsic vorticity of the 
fluid, which is defined as the curl of the generalized canonical momentum and 
which, according to the generalized Larmor Theorem derived in the preceding 
paper, is to be associated with that part of the fluid rotation that is a residual of 
the initial conditions of the fluid. The vorticity invariants are both constants of 
the fluid motion. 

In the case of adiabatic flow, it is possible to express the thermal 4-potential 
in terms of the specific entropy and the temperature integral, which is defined 
as the scalar function whose substantial time derivative is equal to the temper- 
ature. This allows a simple interpretation of the generalized canonical momen- 
tum in terms of the heat reservoir model introduced in the first paper of this 
series. 

iii 
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RELATIVISTIC CHARGED FLUID FLOW 
HI. GENERALIZED HAMILTON-JACOBI EQUATION 

I. INTRODUCTION 

In the preceding paper,* henceforth referred to as II, it was shown that in 
the relativistic dynamics of a charged fluid , a fundamental role is played by the 
intrinsic vorticity tensor where fi is the variable (but relativistically 

invariant) particle mass defined in terms of the constant particle mass m, the 
gravitational potential G , and the specific enthalpy h by the relation 

/jl = m (1 + G/c 2 + h/c 2 ). (1.1) 

From the generalized Larmor Theorem that was derived in II, it was shown that 
the antisymmetric tensor may be regarded as describing the intrinsic rota- 
tion of the fluid, i.e. that part of the total fluid rotation D 1 k that is not produced 
by the action of external fields, but rather is to be associated with the starting 
conditions of the fluid at some instant of past time. 

in 

The physical importance of the intrinsic vorticity tensor 2 /j.co> k liea/the fact 
that the 3-vector given by its space-space components is frozen into the fluid and 
carried along with it. This tensor is defined in terms of the generalized canoni- 
cal momentum by the relation 

2 / no)i k = _ (3i ja k _ 3 k fai), (1.2) 
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where s (j»°> V) is defined by 


jai = /j,yi +(q/c)A j + (m/c)ai; (1.3a) 

or 

c ja 0 = p*c 2 + qA° + ma°, (1.3b) 

= /i*v + (q/c) A + (m/c) a, (1.3c) 

where a j = (a 0 , a ) is a thermal 4-potential that is analogous in its dynamical 
effects to the electromagnetic 4-potential A j s (A 0 , A) , except that the particle 
mass m rather than the particle charge q plays the role of coupling constant. 

As in the preceding paper, the local fluid velocity is designated by vi = T (c, v) 
where c is the speed of light and T = (1 - v 2 /c 2 ) ~ 1/2 where v is the magnitude 
of the 3-velocity v . The mass = T/j. is the mass as seen in the observer's 
frame of reference, as opposed to the invariant mass f± which is the particle 
mass in the fluid rest-frame. 

The equation of motion of the thermal 4-potential is given by (4.9) of II: 

d a Vd t = v k cM a k + (c/m) (ihT^s-tH - 77 i ) (1*4) 

where T and s are respectively the temperature and specific entropy in the 
fluid rest-frame, tt> is the energy -momentum per particle that is injected into 
the fluid because of heat absorption, and 17 > is the viscous 4-force. 

It was shown in the preceding paper that must be orthogonal to v k , 
that is 


a>i k v k = 0. 


(1.5) 
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Contraction of (1.2) with v k yields Euler’s equation: 

d (/xv j )/d r = m B j G + m (B j h - TB j s) + (q/c)F jk v k + 77 j + rj i (1.6) 

where 

F jk = B’ A k - B k A j (1.7) 

is the electromagnetic field tensor. It was pointed out in (5.2) of the first paper 
of this series* (henceforth referred to as I) that 

m(Bi h - TB j s) = p~ 1 3i P (1.8) 

where p is the invariant particle (not mass) density and P is the partial pres- 
sure of the charged fluid. Substitution of (1.8) into (1.6) yields the more usual 
form of Euler's equation. The form given in (1.6), however, has the advantage 
that it does not explicitly involve p, and so has the form of a single-particle 
equation. Once v i has been found by solving (1.6), p can be found from the 
continuity equation: 


d In p/dr = v 1 B. In p = - B. v*. (1.9) 

If we regard the fields G and F ik that appear in (1.6) as given space-time 
functions, then the independent variables of the problem may be taken to be the 
three components of the 3 -velocity v , the particle density p , the enthalpy h, and 
the entropy s - a total of six degrees of freedom. (From the thermodynamic 
properties of the fluid the temperature T is regarded as a known function of the 
variables p, h, and s, and from the physical nature of the heat injection and 
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viscosity tt* and 77 i are regarded as known functions of p, h, s , and v.) The 
corresponding six equations necessary for the solution are the three independent 
space-like components of (1.6) (the time-like component being derivable from 
these), the continuity equation, the thermodynamic equation of state of the fluid 
in terms of p, h, and s, and an equation that specifies ds/dr in terms of the 
six independent variables. 

We could, however, use (1.4) to eliminate s, tt> , and rj 1 from (1.6): 

d (/xv i)/d r = m B * (G + h) + (q/c) F jk v k + (m/c) f jk v fe) (1.10) 

where 

f j k = Bi a k — 3 k a* (1.11) 

is a thermal field that is analogous in its dynamical effects to F> k . This thermal 
field is specified in terms of the 4-potential a i whose "field equations" are given 
by (1.4). Our independent variables could then be taken to be the three components 
of v, the density p, the enthalpy h, and the three significant degrees of freedom 
of ab (One of the four degrees of freedom of a> corresponds to the specifica- 
tion of its gauge, which has no physical significance.) Thus we are now involved 
with a total of eight independent variables (v, p, h, and a). The corresponding 
eight equations are the three independent space-like components of (1.10), the 
continuity equation, the thermodynamic equation of state (which now involves 
a and v as well as P and h), and the three space-like components of (1.4) (the 
time-like component simply specifying the gauge of a j ). 

Comparing these two approaches, we see that the second one, which utilizes 
a, has the disadvantage that it adds two extra variables to the problem. This 
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disadvantage can be removed, however, by taking the fluid equation of motion to 
be a generalized Hamilton- Jacobi equation, rather than the Euler equation. The 
velocity v is then replaced by the scalar Characteristic Function S, so we again 
have a total of only six variables (S, p, h , and a ), even though we continue to use 
a. The desired Hamilton-Jacobi equation is derived in Section II. 

In Section III it is shown that it is always possible to describe the intrinsic 
vorticity 2p.o>i k in terms of two scalars M and 4>, called vorticity invariants, 
that are constants of the fluid motion. The introduction of these two scalars 
leads to a generalized Clebsch Transformation according to which the general- 
ized canonical momentum V is expressed in terms of three scalar functions, 
the Characteristic Function S and the vorticity invariants M and 4>. 

In Section IV the foregoing results are specialized to the case of adiabatic 
(but not necessarily isentropic) flow of an inviscid charged fluid. In this case 
it is possible to express the thermal 4-potential a> in terms of two scalar func- 
tions, one of which is the specific entropy s. The other is a temperature integral 
3 whose substantial time derivative dS/dr is just the local fluid temperature T. 

If c B* 3 is identified with the temperature 4-vector of a heat reservoir 
which coexists in space with the fluid but which, in the adiabatic case, does not 
interact with it, then it can be shown that the entropy-dependent contribution to 
the generalized canonical momentum V may be interpreted in terms of this 
heat reservoir. Because the temperature 4-vector of this reservoir is ex- 
pressible as the gradient of the scalar c3, it is called a scalar reservoir. 
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The general problem of nonadiabatic flow discussed in Section II is formu- 
lated in terms of the six independent variables S, p, h, and a, the last three of 
which are the space-like components of the 4-vector a* . In Section IV, however, 
it is shown that the adiabatic case admits a scalar formulation in that it can be 
formulated in terms of the variables S, p, 3 , s, M, and <t>, all six of which are 
scalar functions. 

In Section V it is shown that, at the expense of increasing the number of in- 
dependent variables, a scalar formulation is also possible even in the most 
general case of nonadiabatic flow. This formulation can be given a natural in- 
terpretation in terms of four imaginary scalar hear reservoirs that exchange 
heat with the fluid. An alternative form of this scalar formulation is also given 
that introduces into the formalism three constants of motion that are related to 
viscosity and heat exchange with the fluid. 

II. GENERALIZED HAMILTON-JACOBI EQUATION 
From (1.2) it is evident that 2/xw jk is expressible as the curl of a 4-vector, 
which we shall designate (m/c) b L It is also evident that this 4-vector must 
differ from -jai by no more than the gradient of some scalar function. Thus 

2 pco> k = (m/c) b k - 3 k b* ) (2.1) 

where 

- (m/c) b^ = jai + S, (2.2) 

where s is an as yet unspecified scalar function. Using (1.3) we have 

- S = pv> + (q/c)A j + (m/c) a > ; (2.3a) 
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or 


- BS/3t = yu.*c 2 + qA° + ma°, (2.3b) 

V S = / i*y + (q/c) A + (m/c) a, (2.3c) 

where 

oi} = (a 0 , a) = a j + b'. (2.4) 

Notation introduced in the preceding two papers has been used in writing (2.3b) 
and (2.3c). In particular 3 j = (3/c Bt, -V); v* = T (c,v) where r = (l _ v 2 /c 2 )” 1/2 
M* s and A j = (A 0 , A) . 

Contracting the curl of a j with v k , and using (1.4), (1.5), and (2.1), we ar- 
rive at the equation of motion of a) : 


da'/d r = v k 3^ a k + (c/m) (m T 3* s - tt> - 17 i ). (2.5) 

This has exactly the same form as the equation of motion of a j given in (1.4). 
From (1.5) and (2.1) it follows that the equation of motion of b* is 

dbi/dr = v k B> b k . (2.6) 

Thus from (2.4) and the forms of the equations of motion of a i, a j , and b j , it 
follows that a* may be regarded as the "driven" part of ai , i.e. the inhomo- 
geneous part of the solution of (2.5), that responds to the "driving force" 

(c/n )(mTB j - 77 j - 7 ?j), whereas b j is the homogeneous part of the solution which 
is insensitive to the driving force. According to ( 2 . 1 ), it is this homogeneous 
part of a) whose curl is the intrinsic vorticity 2/saji k . Thus the description of 
the intrinsic vorticity and the dynamical effects of the viscosity and the entropy 
4- gradient are combined in the single 4-potential ai , which will be called the 
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fluid 4-potential to distinguish it from the thermal 4-potential a 2 , and from b 2 , 
which will be called the vorticity 4-potential . 

In this way we can give a physical significance to the homogeneous part of 
the solution of (2.5), which is not uniquely determined by (2.5) alone, but is also 
dependent upon the initial conditions of the fluid specified at some past instant 
of time. The fact that this homogeneous part of a* determines the intrinsic 
vorticity 2 /xco ik is consistent with the fact, emphasized in the preceding paper, 
that in the generalized Larmor Theorem oo lk is that part of the fluid rotation 
that is determined by the initial conditions of the fluid. This dependence on 
initial conditions will be made very transparent in the next section when it will 
be shown that 2 fj.co* k can be expressed in terms of two scalar functions that are 
constants of the fluid motion. 

It is our objective to arrive at a scalar fluid equation of motion which does 
not involve v 2 , and which has the form of a generalized Hamilton- Jacobi equa- 
tion. Such an equation follows directly from (2.3a) and the velocity normaliza- 
tion condition 

v. v 2 = c 2 . (2.7) 

j 

Taking the terms (q/c) A 2 and (m/c) a 2 to the left side of (2.3a) and equating 
the norms of both sides, we arrive at the desired equation: 

(B. S + qA./c + raa./c) (B 2 S + qA 2 /c + ma 2 /c) = (/xc) 2 . (2.8) 

An alternative form of this equation that makes the effect of a 2 more ex- 
plicit can be derived by taking only (q/c) A 2 to the left side of (2.3a) and equating 
norms of both sides: 
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(3^ S + q A./c) (3 j S + qA j /c) = (/ j . c) 2 + (m/c) 2 a. a j + 2(/am/c)v. a j . (2.9) 

The gauge indeterminacy of a> may be used to simplify this equation by choos- 
ing the gauge so that the condition 

v. =0 (2.10a) 

is satisfied, which implies that 

aj = [(v • a / c), a j. (2.10b) 

This choice of gauge (which is indicated by the overhead tilde) will be called the 
"space-like gauge." Given any a > , the corresponding a } is 

ai = ^ + Bj 0 (2.11a) 

where 0 is the solution of the scalar equation 

d0/dr = - v. aj. (2.11b) 

Note that from (2.10b) it follows that 

a. a j = [(v • a) 2 /c 2 - a • a ] <0. (2.12) 

That is, the norm of aj is always negative. 

When we choose the space-like gauge for aj, the Hamilton- Jacobi equation 


in the form (2.9) becomes 



where the tilde over the S gives explicit recognition to the fact that the Charac- 
teristic Function S which corresponds to a- 1 is in general different from the S 
which corresponds to a' . 

From (2.12) it is evident that the effect of a* is to diminish the right side 
of (2.13), regardless of the nature of the nonisentropy, viscosity, or intrinsic 
vorticity described by at 

As a preliminary to finding the nonrelativistic limit of (2.8), we introduce 
the Principal Function & defined as 

& = S + mV' x. (2-14) 

where x' is the position 4-vector and V' is an arbitrary constant 4-velocity 
which, if the system under study is closed and so conserves its total 4-momentum, 
may be identified with the velocity of the center of mass of the system. From 
(2.3a) and (2.14) it follows that 

- B' & = (/u.v j - mV') + (q/c) A' + (m/c) a.1 (2.15) 

which, together with the normalization condition V. V' = c 2 , yields the following 
form of the Hamilton-Jacobi equation: 

V* 3 j& - (B. & + q Aj /c + ma./c) (3' £> + qA'/c + mai/c) /2 m 

(2.16) 

+ qAj V'/c + ma. V j /c + m(G + h) [1 + (G + h)/2c 2 ] = 0. 

Now we note that in the special frame of reference for which 

V j = (c, 0, 0, 0) (2.17) 
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(2.16) may be written as follows: 


9 1 &/B t + (V^> + q A/c + mo/c) • (V$ + q A/c + ma/c) /2 m 

(2.18) 

+ qA° + ma° + m(G + h) = [(3 &/'d t + qA° + ma°) - (mG + mh) 2 ] /2mc 2 . 

From (1.1) and the time-like component of (2.15) for the frame in which V j has 
the form (2.17), we find 

[(3<?>/5 t + q A 0 + m a 0 ) - (mG + mh) 2 ]/2 m c 2 

= c 4 [(ju* - m) 2 - (/x - m) 2 ]/2mc 2 (2.19) 

£ (1/2) mv 2 [(1/2) mv 2 + 2 (mG + mh) ]/2mc 2 . 

Because this is obviously zero in the nonrelativistic limit, the right side of (2.18) 
vanishes in this limit, and we are left with an equation having the form of the 
familiar nonrelativistic single-particle Hamilton-Jacobi equation, except for the 
additional terms involving h, a 0 , and a, which play roles analogous to G, A 0 , and 
A respectively. 

IH. GENERALIZED CLEBSCH TRANSFORMATION 
It was remarked in the preceding section that the vorticity 4-vector b^ must 
be regarded as a function of the initial conditions of the fluid, and that this could 
be made most evident by showing that b* can be expressed, quite generally, in 
terms of two scalar functions that are constants of the fluid motion. Thus the 
specification of these scalars everywhere at some instant of past time suffices 
to determine them, and hence b j , for all time. 


11 



As a preliminary to proving this statement, we first recall that in (4.12) of 
II it was pointed out that the orthogonality condition co ,k v k = 0 requires that co ik 
have the following form: 

(a) 10 , a; 20 , o> 30 ) = - r«o X v / c, (3.1a) 

(co 23 , w 31 , w 12 ) = r «, (3.1b) 

where (3.1b) constitutes the definition of the 3-vector «. It was further pointed 
out in (4.18) of II that a necessary consequence of the requirement that co> k have 
the form given in (3.1), and the requirement that 2/u.ad k be expressible as the 
curl of a 4-vector (cf. (2.1) above), is that the 3-vector intrinsic vorticity 2/x*« 
must satisfy the following two equations: 

V • (2/la*co) = 0; (3.2a) 

B(2/x*«)/3t = Vx [vx (2/n*»)]. (3.2b) 

For the moment let us confine ourselves to a fixed frame of reference, and let 
us assume that 2 is continuously differentiable. It is well known 1 that such 
a vector, which by (3.2a) must be solenoidal, can always be expressed as follows 
in terms of two scalar (in the three-dimensional sense) functions, which we shall 
denote M and 4> : 


2/^.*w = (VM) x (V0). (3.3) 

For given 2 /**<*, this equation which follows from (3.2a), may be regarded as a 
condition on the spatial dependence of M and 4>. Eq. (3.2b) represents a condi- 
tion on both the spatial and time dependence of M and <f>. This condition is equiv- 
alent to the following requirement: 
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v x (2 / 2 *w) = (3 M/3 t) - (3<£/3 t) VM. 


(3.4) 


That (3.4) is equivalent to (3.2b) is easily verified by taking the curl of (3.4) and 
using (3.3). The gradient of an arbitrary scalar function could have been added 
to the right side of (3.4) without affecting its equivalence to (3.2b), but because 
this scalar function is arbitrary, we are free to equate it to zero with no loss of 
generality. 

We have shown that, because 2 m u > jk is the curl of a 4-vector, 2 w must be 
expressible in terms of two scalars M and 4 > that satisfy the condition (3.4). So 
far we have restricted ourselves to a single frame of reference, and have said 
nothing about the Lorentz transformation properties of 2 l tx*w. The fact that this 
3 -vector is the space- space part of an antisymmetric world tensor represents an 
additional restriction on the space-time dependence of 2 m*“, and hence of M and 
4 >. This requirement on M and 4 > turns out to be just that they both be scalars 
not only in the three-dimensional sense, but also in the four-dimensional sense 
as well. That is, they must be Lorentz invariant. The validity of this statement 
follows from the fact that, making use of (3.1), we can write the six component 
equations of (3.3) and (3.4) in the following way, which has the formal appearance 
of a tensor equation: 


2 A^ jk = (3 j M) (3 k 0) - (3 k M) (3 j 4 >) 
= 3 j (M 3 k 0) - 3 k (M3 j 4 >). 


(3.5) 


If now we require that M and 4 > be Lorentz invariant, then the right side of (3.5) 
is a genuine world tensor, and consequently so is 2fto> jk Q n the left side, which 
guarantees that 2^*« has the correct Lorentz transformation properties. 
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There is a more abstract, but much more direct, way of arriving at (3.5). 
Referring to (4.18) of II, we note that the four component equations of (3.2) can 
be written in the following form: 

d> (2 AioX) + 3 k (2 /J-ccfi j ) + 3^(2Ai.oi k ) = 0. (3.6) 

The tensor 2/i^ ik has six degrees of freedom (before we impose the orthogonality 
requirement (1.5)). But (3.6), which is equivalent to the four component equations 
of (3.2), removes four of these degrees of freedom. Thus any antisymmetric 
world-tensor with two degrees of freedom that automatically satisfies (3.6) is a 
perfectly acceptable and general way of representing 2^u lV . It is obvious that 
(3.5) satisfies these requirements, which is all the justification it needs. 

Having derived (3.5) from the requirement that 2(UO) jk be the curl of a 4- 
vector, we may now impose the orthogonality condition (1.5) in order to derive 
the equations of motion of M and <£. Contracting (3.5) with v k and using (1.5) 
we find 

(3 > 4>) d M/d r = (3 i M) d <p / d r. (3.7) 

Constructing the four-dimensional cross-product of this equation with 3 k M and 
using (3.5), we have 


2/xa)i k dM/dr = 0. (3.8) 

If 2/u.o>i k = o, then from (3.5) it is evident that M and 4> may be set equal to con- 
stants, and so may be considered constants of the motion. In the non- trivial case 
in which 2 /ua>i k ^ o , it follows from (3.8) that 

dM/dr = 0, (3.9) 
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which says that M is a constant of motion. Similarly, by constructing the cross- 
product of (3.7) with B k 0, we arrive at the conclusion that 

d0/d r = 0. (3.10) 

Thus both M and 0 are constants of motion. 

Comparing (2.1) with (3.5), we arrive at the following expression for the 
vorticity 4-vector b j in terms of the vorticity invariants M and 0 : 

(m/c)^ = M3 j 0, (3.11) 

where we note that the gradient of an arbitrary scalar function might have been 

added to the right side of (3.11), but because only the curl of bi is of physical 

significance, this scalar function has been set equal to zero, with no loss of 

generality. It is easily verified that the equations (3.9)-(3.11) guarantee that the 

Sst'.S-fl frd, 

equation of motion for bi given in (2.6) is automatically ynrtfiwd . 

Substituting (3.11) into (2.2), we arrive at the generalized Clebsch 
Transformation: 

P = - (B j S + M3 j 0); (3.12a) 

or 

cfa° a _ (3S/Bt + MB0/3t), (3.12b) 

= VS + MV0, (3.12c) 

where the definition of the generalized canonical momentum & = (fc°. fc) is 
given by (1.3). Unlike the familiar form 2 of the Clebsch Transformation, which 
expresses the fluid velocity in terms of three scalars, the generalized trans- 
formation refers to the generalized canonical momentum. 
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An important special case arises if 2p.*« = 0 . From (3.2b) it is evident that, 
if this condition holds everywhere in 3 -space at any instant of time, then it will 
propagate itself for all time. From (3.3) or (3.5), it is evident that in terms of 
M and 4 > this zero-vorticity case corresponds to one of the following conditions: 
Either one or both of the scalars M and 4 > is everywhere constant for all time, 
or else M can be expressed as a function of 4 >. When any of these conditions 
holds, it is evident that the term -MB’ 4 > in (3.12a) can be expressed as the gra- 
dient of a scalar function, which function could then be absorbed into the func- 
tion S. Thus the zero-vorticity case can be characterized by the condition 

(Potential Flow) = - B j S, (3.13) 

and so will be referred to as generalized potential flow . It is the self-perpetu- 
ating feature of this special type of flow that makes it important. This self- 
perpetuation is intuitively obvious in terms of M and <£: Because these are con- 
stant along every particle trajectory, if either of them is the same for all 
trajectories at any instant of time, then it must be the same for all time. The 
formal expression of this remark for the case of M, for example, follows from 
(3.9) which can be written 


BM/Bt = -,v • VM- (3.14) 

It is obvious that if VM= 0 , then BM/B t = 0 and the spatial uniformity of M 
perpetuates itself for all time. 

Our intuitive feeling for the physical meaning of M and 4 > can be helped by 
noting that, because w is an angular velocity, it follows from (3.3) that the 
product must have the dimensions of angular momentum. With no loss of 
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generality we may assert that <p is dimensionless and that M has the dimensions 
of angular momentum. Thus M may be regarded as an angular momentum (per 
particle) - in some way related to the intrinsic vorticity - which is conserved 
along the particle trajectory. We may think of <p as the initial value at t = 0 of 
one of the particle coordinates (in dimensionless form), whose memory is re- 
tained by the particle for all time. 

A simple example that illustrates these points is provided by the case of 
nonrelativistic rigid rotation of the fluid in the absence of any external fields. 

The rotation is thus characterized by a single angular velocity vector of mag- 
nitude fi that is everywhere the same. It is easily shown 3 that in this case M 
turns out to be equal to the angular momentum of a particle about the rotation 
axis, i.e. M = mfir 2 where r is the particle distance from the axis. 4> turns out 
to be the value cp 0 of the particle’s azimuthal coordinate at t = 0 , i.e. </> = cp-ft t, 
where the azimuthal angle <jp is given by cp = ft t + cp Q . 

If there were a uniform magnetic field parallel to the rotation axis, then the 
expression for <j> would remain unchanged, but M would become M = m(fi-Q L ) r 2 - cor 2 
where fi L is the Larmor rotation velocity produced by the magnetic field, and 
co = Cl - 0 L is just the magnitude of the intrinsic angular velocity. Thus M would 
be that part of the total particle angular momentum that is produced by the in- 
trinsic angular velocity co. 

Note that in these examples the existence of a nonzero intrinsic angular 
velocity is associated with an M that is not spatially uniform. Tills is in fact 
always the case because, if VM = 0 , then from (3.3) it follows that « = 0 and, 
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according to (3.14), once V M = 0 , this condition will perpetuate itself. It is in- 
teresting to speculate 3 that there might exist a natural mechanism, namely 
turbulence, that would tend to make M uniform even if this were not originally 
the case. The argument goes as follows: Inasmuch as M is a constant of mo- 
tion, it may be regarded as an intrinsic property of every small sample of the 
fluid, which is carried along with the fluid and characterizes it for all time. 
Because turbulence tends to mix the fluid, it is reasonable to expect that it might 
produce a diffusion of the nonuniformities in M , with the result that M would tend 
to become uniform throughout the fluid (to the extent that it is not impeded from 
doing this by imposed constraints and boundary conditions 3 ). 

Because uniform M implies the existence of generalized potential flow as 
defined by (3.13), the above speculation suggests that potential flow may be en- 
countered very frequently in situations characterized by strong turbulence, i.e. 
by a very large Reynolds number. 

IV. ADIABATIC FLOW 

In the preceding section it has been shown that it is always possible to ex- 
press the vorticity 4-potential b j in terms of two scalar functions. It will now 
be shown that, in the case of adiabatic flow, it is possible to express the thermal 
4-potential a j in terms of two scalar functions, one of which is the specific 
entropy s . 

Adiabatic flow is characterized by the fulfillment of the following three 
conditions: 

§ s d s/dr = 0; (4.1a) 
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7r J = 0; 


(4.1b) 


rji = 0. (4.1c) 

From the physical point of view these conditions are not independent, because 
the first implies the other two. Thus, if its local entropy cannot change, the 
fluid cannot (reversibly) gain or lose heat, which implies that tt' } = 0 . As to 
(4.1c), which says that the viscous force must be zero, this follows from the ob- 
servation that viscosity always causes entropy generation. 

From the formal point of view, however, if we are not given the functional 
dependence of tt 1 and 17 ) on s, the conditions (4.1b) and (4.1c) must be regarded 
as independent of (4.1a). Such a functional relation for tt> was in fact given in 
(5.7) of I: 


7 t> = mTjs/c (4.2) 

where is the temperature 4-vector of the heat reservoir from which heat is 
reversibly transferred to the fluid. The formal condition guaranteeing reversi- 
bility of the transfer was given in (4.27a) of I: 

T = v. T jj/ c. (4.3) 

(Following the notation change noted in II, we designate the fluid temperature in 

0 

the fluid rest-frame by T rather than by T as in I.) It is obvious that, given 
(4.2), the condition (4.1b) follows from (4.1a). Lacking an expression for analogous 

to the one for tt* given in (4.2), we must regard (4.1c) as being formally independ- 
ent of (4.1a). In any case, we assume that all three conditions of (4.1) are fulfilled. 
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Thus (1.4) may be written in the form 


(B j a k - B k a j ) v k = - cT3's. (4.4) 

The solution of this equation can be expressed directly in terms of s and the 
temperature integral 4 0 which is defined by the relation 

T s d3/dr s v. B j 3. (4.5) 

It is easily verified that when (4.1a) holds, either of the following expressions 
for a> satisfies (4.4): 


a j = - c s B 1 3; 
a* = c3B> s; 

which proves our assertion that for adiabatic flow it is possible to express 
in terms of two scalar functions, one of which is s. 


(4.6) 

(4.7) 


a* 


Note that, because these two expressions for the thermal 4-potential differ 
only by the gradient of c sO, they are actually the same 4-potential with different 
choices for the gauge. Because of (4.1a), a j is obviously orthogonal to v. : 

v. a ‘ = c 3 s =0. (4.8) 

j 

Although this orthogonality is a formal advantage, the choice of gauge made in 
writing (4.6) has the intuitive advantage that it lends itself to an explanation of 
the physical significance of a j , or more particularly, of the term (tn/c) a> that 
appears in the expression for the generalized canonical momentum that is 
given in (1.3a). 
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Intuitive Interpretation of Generalized Canonical Momentum 

This intuitive interpretation is based on an identification of c with the 
heat reservoir temperature 4-vector T^, which is suggested by a comparison 
of (4.5) with (4.3). Thus we assert that is defined by the relation 

(4.9a) 

R 

Because T ^ = T° (1, v R /c) where T°, the time-like component of T£, is the 
reservoir temperature as seen in the observer's frame, and v R is the reservoir 
3 -velocity, it is possible to write (4.9a) in the following form: 

= (4.9b) 

T rV c2h - V3 - (4,9c) 

Whether or not the heat reservoir specified by may actually be identified 
with whatever physically real reservoir that may happen to coexist in space with 
the fluid (but not interact with it for the case of adiabatic flow) depends on whether 
or not the temperature 4-vector of this real reservoir can be expressed as the 
gradient of a scalar. (Such a reservoir will henceforth be called a scalar reservoir .) 
A restriction of this kind can be expected to exclude many physically realizable 
reservoirs. Thus the intuitive argument given below, which is based on the identi- 
fication made in (4.9), has physical significance only for this rather narrow class 
of scalar reservoirs. 
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Using (1.1), (4.6), and (4.9) in (1.3), we find 
= ^v 1 + (q/c)A j -(m/c)sT^ 

= mv j + (mG/c 2 ) v j + (q/c) A j + m [ (h/c 2 ) v j - s T^/c]; 
or 

c)a 0 = p.*c 2 + qA° - m s T R 

= m*c 2 + m*G + qA° + (m*h - m s T R ), 


(4.10a) 


(4.10b) 


^ - fj.*\ + (q/c) A - (m sT“/c 2 ) v r 

(4.10c) 

= m*v + (m*G/c 2 ) v + (q/c) A + [ (m*h/c 2 ) v - (ms T R /c 2 ) v R ] . 


Note that by contracting (4.10a) with v., we have 

c'ta°=mc 2 +mG + qA°+m(h-sT), (4.H) 

O O 

where p° = v. f> J /c and A 0 = v. AVc are respectively the time-like components 
of and A j as seen in the fluid rest-frame. 

Thermodynamics is represented in (4.11) by the specific Gibbs function 
g = h - sT- From the intuitive point of view, this is not a surprising result. 
Each small sample of fluid, because it is constrained to have the same pressure 
and temperature as the surrounding fluid, seeks a thermodynamic equilibrium 
characterized by the constraint that virtual displacements from equilibrium 
must produce no change in either pressure or temperature. But the appropriate 
thermodynamic potential function in such a case is the Gibbs function. Thus it 
is not surprising that in (4.11), which gives the total rest-energy of a small 
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sample of the fluid (on a per-particle basis), the thermodynamic energy should 
be represented by the Gibbs function. 

It is possible to use the heat reservoir concept to give an alternative ex- 
planation of the fact that it is the Gibbs function that plays the role of thermo- 
dynamic potential. For simplicity, we shall first consider the case for which 
v R = v , and shall work in the common rest-frame of the fluid and reservoir. 

We first note that the specific enthalpy h, often called the heat content of 
the gas, is just the quantity of heat that would have to be injected into a sample 
of unit mass of the fluid in order to expand it against the pressure P of the sur- 
rounding fluid (which, assuming the sample to be very small, would remain con- 
stant during the expansion) while bringing its temperature from absolute zero up 
to the temperature T of the surrounding fluid. This is obvious from the relation 

o o 

h = u + P l) where u and l) are respectively the internal energy and volume of 
unit mass of the fluid. Of the total heat energy h injected into the sample of 

O 

fluid, the amount PI) is converted into the mechanical work necessary to push 
back the surrounding fluid as the sample of unit mass is expanded from zero 

O 

volume to its final volume l), and the amount u remains in the sample as its 
internal energy. The energy u obviously resides within the fluid sample. The 

O 

work PI), however, was performed on the fluid surrounding the sample, and so 
is a potential energy that is stored outside of the sample. Because, however, 
the heat energy h was injected into the space occupied by the sample in expand- 

O 

ing it to its final volume i), ana could in principle be extracted from this region 
of space by reversing the process, it is legitimate to associate the energy h 
with the region of space occupied by the fluid sample. But now we recall that. 
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from the point of view of the heat reservoir model, this region of space is occu- 
pied by, both the fluid and the reservoir , which are regarded as separate systems. 
The question thus arises as to how much of the energy h should be associated 
with the fluid, and how much should be assigned to the reservoir. It is easy to 
see that this latter energy is just s T (per unit mass of the fluid). The reason 
for this is that, if an external reservoir at absolute zero of temperature were 
available, a Carnot engine could extract the energy s T from the fluid reservoir. 
Thus the energy to be associated with the fluid is h - sT, which is just the Gibbs 
function g. It is to be expected that, in a dynamical theory concentrating on the 
fluid , rather than on the heat reservoir associated with it, the thermodynamical 
energy is represented by the fluid energy g = h - s T rather than by the total 
energy h. The relation (4.11) confirms this expectation. 

This argument is easily generalized to the case for which v R ^ v . In this 
case, in the fluid rest- frame we may still associate the energy h with the region 
occupied by unit rest-mass of the fluid. In the observer's frame, however, this 
energy becomes Th, and necessarily associated with it is the momentum (rh/c 2 )v. 
Thus we associate the energy-momentum 4-vector (h/c 2 )vj = [(rh/c), (fh/c 2 ) v] 
with the volume occupied by unit rest-mass of the fluid. 

Next, referring to (4.10) of I, we note that, if an external reservoir at abso- 
lute zero of temperature were available, a Carnot engine could extract from the 
fluid reservoir (per unit rest-mass of fluid) an amount of energy, and momentum 
given by the 4-vector (s T ^/c) = [(s T R /c), (s T°/c 2 ) v R ] , where we have used 
the fact that = T°(l, v R /c ). Thus the energy-momentum to be associated 
with the fluid above, on a per-particle basis^is m [ (h/c 2 ) v j - s T^/cJ. This is 
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just the expression that in (4.10a) represents the thermodynamical contribution 
to the generalized canonical momentum 1» J . 

Finally, it should be noted that we could have inserted a * = c 3 3 J s rather 
than a i - -c sB j 3 into (1.3), in which case, instead of (4.10a), we would have 
arrived at the following expression for the generalized canonical momentum 5 : 

= fx v j + (q/c) A j +m3 3 j s. (4.12) 

Because a 1 and a3 differ only by a choice of gauge, which can have no physical 
significance, whether we use or t in the overall formalism can lead to no 
physically observable differences. Because, as indicated in (4.8), is orthog- 

*v> : . 

onal to v , it turns out that ja has certain formal advantages over V . The 
j 

latter, however, has the intuitive advantage that it lends itself to the simple 
physical interpretation discussed above. 

Hamilton-Jacobi Equation 

The generalized Hamilton-Jacobi equation was given in (2.8) in terms of the 
fluid 4-potential cd and the corresponding Characteristic Function S, and also in 
(2.13) in terms of a. 1 and the corresponding S. The 4-vectors ad and 33 dif- 
fered only in choice of gauge. In the case of S3, the gauge (the space-like gauge) 
was chosen so that the orthogonality condition v. a j =0 was fulfilled. As stated 
in (2.12), a necessary consequence of this choice is the fact that the norm of ai 
is always negative, i.e. a. a' < 0. In the case of adiabatic flow, it is evident 
from (2.4), (4.6), (4.7), and (3.11) that cd and 33 may be written as follows: 

(m/c) cd = -ms3^3 + MJ^</>; (4.13a) 
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(m/c) “a) - m 0 3* s - 4> B* M. 


(4.13b) 


Obviously, these two forms of the fluid 4-potential differ only by the gradient of 
the scalar function (ms 3 - M 4>). Because of (3.9) and (4.1a) it follows that a' 
satisfies the orthogonality condition v. a j = 0 . The form of the fluid 4-potential 
given in (4.13a) is the sum of the expressions for a j and b’ that we have used 
for intuitive purposes. The form a> given in (4.13b), however, has certain 
formal advantages arising from the orthogonality of a* and v.. Moreover, it 
has the advantage that it exhibits more explicitly than ai the fact that if s and 
M become everywhere constant, i.e. if d’ s = M = 0, then the fluid 4-potential 
drops out of the picture. (For this case ai vanishes, whereas ai becomes the 
gradient of a scalar.) For these reasons we shall restrict ourselves to a j and 
the corresponding Characteristic Function S. 

Using (4.13b) in (2.13), we see that the Hamilton- Jacobi equation for adia- 
batic flow becomes 

(B j S + qA j / c) (9 j S + qA j /c) = (/lc c ) 2 + (m 3 3 s - 4> M) (m 3 s - 0^M). (4.14) 

Because of the fact that a .33 < 0, the effect of nonuniformity in s and M is 
always to diminish the right side of (4.14). 

We now augment (4.14) with the equations of motion of 3, s,4>, and M, and 
the continuity equation (1.9), which is the equation of motion of p: 


d3/d r = T; 

(4.15) 

d s/dr = 0 ; 

(4.16) 
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dcfr/dr - 0 ; 


(4.17) 


dM/dr = 0; 

(4.18) 

d p/dr - - p 3. v* . 

(4.19) 


The 4-velocity vi which appears explicitly in (4.19) and implicitly in the time 
differentiation operator d/d r = v j 3. is given by (2.3a): 

V* = -/X' 1 [Bis + (q/c) A j + m3 B> s - ^ B> M] . (4.20) 

We regard G (which appears in /z ) and A j as given space-time functions. If the 
thermodynamic equation of state of the fluid is given, then h (which appears in /z ) 
and T (which appears in (4.15)) are known functions of p and s. (Knowing this 
functional dependence of T, we may regard (4.15) as the fluid equation of state in 
terms of 3.) Thus we have expressed the problem in terms of the six scalars 
S,3, s, M> <£> and p which are determined by the six scalar equations (4.14)-(4.19). 
This could be called the scalar formulation of the dynamical problem to distin- 
guish it from the usual formulation based on the 4-vector Euler equation. 

The remark that, once the fluid equation of state is known, h and T may be 
regarded as known functions of p and s can be illustrated by the case of a per- 
fect gas. In this case 

( h = c p T, (4.21a) 

(Perfect Gas) < where 

t T = K p y ' 1 e rs/c P, (4.21b) 

where K is an arbitrary constant, c p is the constant-pressure specific heat, and 
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y - c p /c v is the ratio of constant-pressure and constant- volume specific heats. 
Using (4.21b) and (4.19), it is easily shown that for adiabatic flow of aperfectgas 
(4.15) may be replaced by the following equation of motion for 3 , which may also 
be regarded as the fluid equation of state in terms of 3: 

(Perfect Gas) d 2 3/dr 2 + [ (y - 1) 3. v j ] d 3/d r = 0. (4.22) 


Three important special cases should be noted. The first of these is 
isentropic flow which arises if 3. s = 0. In this case s drops completely out of 
the formalism. The second case is potential flow which, as noted in Section III, 
arises if either M or 4> becomes constant, or if either can be expressed as a 
function of the other. In this case both M and <f> may be dropped from the for- 
malism. The third special case, isentropic potential flow is simply the simul- 
taneous fulfillment of the above two conditions. In this case the six unknowns of 
the problem reduce to the three unknowns S, 3 , and p. Because 3 appears only 

r\j 

in h, it can be replaced by h , so that the three unknowns could be taken to be S , 
h, and p . For a perfect gas , the system of equations for isentropic potential flow 
becomes 


(c 3. S + qA.)(c 3 j S + qA j ) = (me 2 + mG + mh) 2 ; (4.23a) 


( I sent ropi c 
Potential <( 
Flow) I 


dp x dr = - p Bj v^ ; 


(Per f ect gas) dh/ dr = - (y - 1) h v j ; 
where 


vVc = - (c S + qA j )/(mc 2 + mG + mh). 


(4.23b) 

(4.23c) 


(4.23d) 
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Finally, it should be noted that, although in the system of equations (4.14)- 
(4.19) we regard p as one of the independent variables, it is in fact possible to 

r\j 

regard only the five variables S, 3 , s ,<p, and M as constituting the complete set 
of independent variables. From this point of view, the two independent thermo- 
dynamic variables are taken to be 3 and s, and then h and T are regarded as known 
functions of these (assuming that the thermodynamic properties of the fluid are 
completely known). The thermodynamic equation of state in terms of 3 is given 
by (4.15), where T is regarded as a known function of 3 and s. This equation 
would then involve the thermodynamic variables 3 and s, but not p . (The fact 
that this equation of state involves only two thermodynamic variables is a con- 
sequence of the fact that the flow is constrained to be adiabatic.) From (4.22), 
we see in fact that for a perfect gas it does not even involve s (except for the 
implicit dependence involved in v J that is given by (4.20)). Once the problem 
has been solved and S, 3, s, 0, and M are known space-time functions, p can be 
found from (4.19). 

V. FOUR-RESERVOIR MODEL OF NONADIABATIC FLOW 
In the case of nonadiabatic flow, all entropy-dependent effects (including the 
effects of tt i and rfi which are always associated with entropy generation) are 
represented by the thermal 4-potential a j whose dynamical role is analogous to 
that of the electromagnetic 4-potential. The "field equation" for a j was given 
in (1.4), and may be written 

(m/c) (3 J a k - a j )v k + mT 3 J s = 7 r J + 77 * (5.1) 
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where from (4.2) and (4.3) 


Vj 7r j = m T s (5.2a) 

and from (4.3) of II and (5.4) of I 

v. r)i = 0. (5.2b) 

The 4-vectors n> and if are regarded as known functions of p, , s, and T 
that satisfy the conditions (5.2). The component equations of (5.1) are the four 
equations whose solution yields the space-time dependence of the four compo- 
nents of a’ . 

In addition to this 4-vector treatment of the entropy-dependent forces in the 
case of nonadiabatic flow, it would be desirable to have a scalar formulation that 
was a generalization of the one given in the preceding section for the case of 
adiabatic flow. We arrive at such a generalization by assigning to a'i the func- 
tional form 

BJ 3, (5.3) 

(N) V ' 

(5.4a) 

^(i> -^(2) = ^(3) =3(4) (5.4b) 

where 3 (N) = d3 /dr. The form assigned to a i in (5.3) is an obvious gener- 
alization of the form given in (4.6) in the case of adiabatic flow. The conditions 
given in (5.4) represent four constraints on the eight variables s (N) and 3 (N) 


where 


-z: 


(N) 


Nsl 


* 

L 


S (N) S 
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so the right side of (5.3) has four degrees of freedom which suffice to express any 
arbitrarily given a j . We need four more equations in addition to (5.4) in order to 
arrive at a complete specification of the eight variables. The additional four equa- 
tions are supplied by the following 4-vector equation: 

m V"’ S (N) B j 3 (n) = -rr> + 17 i. (5.5) 

In order to demonstrate that fulfillment of the conditions (5.4) and (5.5) does in- 
deed guarantee that (5.3) satisfies (5.1), we first note that by contracting (5.5) 
with v. and using (5.2) and (5.4), we arrive at the following result: 

3 (n) s d3 (N) /dr = T; N = 1, 2, 3, 4. (5.6) 

It is now easily verified that when (5.4) and (5.5) (which implies (5.6)) are satis- 
fied, then (5.3) does indeed satisfy (5.1). Thus we have replaced the four variables 
ai and the corresponding four equations (5.1) with the eight scalar functions s (N) 
and c^ N) and the corresponding eight equations (5.4) and (5.5). 

In terms of number of variables, this replacement is obviously not advan- 
tageous. The fact that the new variables are scalars rather than the components 
of a 4-vector can, however, be an important advantage. Moreover, it is possible 
to give the scalar formulation a simple physical interpretation in terms of. a 
thermal interaction of the fluid with four separate heat reservoirs. The basis of 
this interpretation is the fact that, because of (5.6), it is possible to identify each 
of the four gradients cS j 3 with the temperature 4-vector T j N) of a scalar heat 
reservoir in thermal interaction with the fluid: 


Bi 3 s Ti 

w w / xi \ 1 / 


<N) 


(N)' 


(5.7) 
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This is an obvious generalization of (4.9a). The 4-vector equations (5.3) and 
(5.5) may now be written as follows: 


(m/c) a } 


s T 1 
(N) A (N)’ 


(m/c)^ s (N) T} N) =tt i + V >. (5.9) 

There is a simple physical interpretation of (5.9). Referring to (4.22) of I ; we 
note that m s (N) T J (N) may be interpreted as the energy-momentum (per particle) 
delivered by the Nth heat reservoir per unit time to the fluid, and is that 

part of s, the total fluid entropy increase per unit time, that results from the 
absorption of heat from the Nth reservoir. Note that it is even possible to ac- 
count for a purely space-like viscous force on the right side of (5.9) (i.e. = 0) 

in terms of this four-reservoir model. In such a case, for example, we might 
have two reservoirs moving in opposite directions, one absorbing and the other 
rejecting heat in such a way that no net heat energy is delivered to the fluid (in 
its own rest-frame), but a net momentum per unit time (i.e. force) would be 
delivered. This force would be the viscous force. In such a case we would have 
s (1) > 0 , s (2) < 0, and s = s (1) + s (2) = 0 , so the fluid flow would actually 
be adiabatic even though a viscous force was present. (Adiabatic flow of a viscid 
fluid is, of course, a mathematical idealization in that physically viscosity is 
always accompanied by entropy generation.) 


Substituting (5.8) into (1.3a), we arrive at the following expression for the 
generalized canonical momentum in the case of nonadiabatic flow: 
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(5.10) 


= + (q/c) A J 


-*z 


5 T 1 
(N) A (N) 


Referring to the interpretive discussion given in the preceding section, we note 
that ms (N) may be interpreted as the energy-momentum (per particle) 
contained in the Nth reservoir, which could in principle be converted into 
mechanical form by a Carnot engine operating between the Nth reservoir and 
a cold reservoir at absolute zero of temperature. By subtracting from the total 
energy-momentum /^v^ + (q/c) A j that part ; m 2 s (N) T which should be asso- 
ciated with the four reservoirs, we are left with the part ft , which should be 
associated with the fluid. 


Although the above formulation of the nonadiabatic flow problem lends itself 
most readily to physical interpretation, there exists an alternative formulation 
that has certain formal advantages. We arrive at this alternative by writing 
as 

5 (n) = + ^ N) (5.11a) 

where 

3' ■ dS'/dr = T (5.11b) 

and 

F' N) = dF' N) / dr = 0, (5.11c) 


where N = 1, 2, 3, 4. The fulfillment of these conditions will automatically guar- 
antee the fulfillment of (5.4b) and (5.6). Substituting (5.11) into (5.3) and using 
(5.4a), we may write 


a j 


c 


s a* 3 + L S (N) ^ F (N) 


(5.12) 
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where 


3 = 3 ‘ +«"(«>; F ( „ =f;„,-fj 4) , (n = i, 2, 3, )- 


The conditions on 3 and F (N) that correspond to (5.11) are 


3 = T; 


F = F = F =0 

O) " (2) r (3) U ’ 


(5.13) 


(5.14a) 

(5.14b) 


These are four conditions on the eight variables s, 3 , F (n) , and s (N) (N = 1,2,3). 
The remaining necessary four equations are provided by the following 4-vector 
equation: 


m s (N) F (n) = V +r}> -ms3i 3 . (5.15) 

N ” 1 

It is easily verified that if (5.14) and (5.15) are satisfied then (5.12) satisfies 
(5.1). 


The condition (5.14a) specifies 3 only to within an arbitrary additive func- 
tion that is a constant of motion. We could use this freedom to specify 3 so that 
m sB^3 is the closest possible approximation to rr j + t]\ Then we would interpret 
the left side of (5.15) as expressing that part of tt* + 77 j that is not describable 
in terms of a single scalar reservoir model. The formal advantage of (5.12) 
over (5.3) is that -n-i + 17 j is described in terms of three constants of motion 
F (N) and corresponding entropy contributions s (N) (N = 1,2,3) upon which no 
constraint of the kind given in (5.4a) is imposed. Moreover, the fluid equation 
of state can be written simply in terms of p, 3, and s alone, making it unneces- 
sary to use all four entropy contributions s , x as would be the case for the 

(N) 

formulation based on (5.3). 
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VI. CONCLUSIONS 


Use of the generalized Hamilton- Jacobi equation for a charged fluid in the 
presence of given electromagnetic and gravitational fields allows the problem 
to be formulated in terms of a set of independent dynamical and thermodynamical 
variables that are all scalars. Such a scalar formulation is possible even when 
viscosity and heat injection into the fluid are taken into account. 
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